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ABSTRACT 

Line-driven winds from hot stars and accretion disks are thought to adopt a unique, critical 
solution which corresponds to maximum mass loss rate and a particular velocity law. We show 
that in the presence of negative velocity gradients, radiative-acoustic (Abbott) waves can drive 
shallow wind solutions towards larger velocities and mass loss rates. Perturbations introduced 
downstream from the wind critical point lead to convergence towards the critical solution. By 
contrast, low-lying perturbations cause evolution towards a mass-overloaded solution, developing 
a broad deceleration region in the wind. Such a wind differs fundamentally from the critical 
solution. For sufficiently deep-seated perturbations, overloaded solutions become time-dependent 
and develop shocks and shells. 

Subject headings: accretion disks — galaxies: active — cataclysmic variables — hydrodynamics 
— stars: mass loss — stars: winds 
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1. Introduction 

Atmospheres of hot luminous stars and accretion 
disks in active galactic nuclei and cataclysmic vari- 
ables form extensive outflows due to super-Eddington 
radiation fluxes in UV resonance and subordinate 
lines. An understanding of these winds is hampered 
by the pathological dependence of the driving force 
on the flow velocity gradient. Castor, Abbott, & 
Klein (1975; CAK hereafter) found that line-driven 
winds (hereafter LDWs) from O stars should adopt 
a unique, critical state which corresponds to maxi- 
mum mass loss rate. The equation of motion for a 1- 
D, spherically symmetric, polytropic outflow subject 
to a Sobolev line force allows for two infinite fam- 
ilies of so-called shallow and steep solutions. How- 
ever, none of these families can provide for a global 
solution alone. Shallow solutions do not reach in- 
finity, while steep solutions do not extend into the 
subsonic regime including the photosphere. The crit- 
ical wind starts then as the fastest shallow solution 
and switches at the critical point in a continuous and 
differentiable manner to the slowest steep solution. 
Hence the critical point and not the sonic point de- 
termines the bottleneck in the wind. This description 
in principle applies equally to winds from stars and 
accretion disks. 

A physical interpretation of the CAK critical 
point was given by Abbott (1980), who derived a 
new type of radiative-acoustic waves (hereafter Ab- 
bott waves). These waves can propagate inward, in 
the stellar rest frame, only from below the CAK criti- 
cal point. Above the critical point, they are advected 
outwards. Hence, the CAK critical point serves as 
an information barrier, much as the sonic or Alfven 
points in thermal and hydromagnetic winds. Ab- 
bott's analysis was challenged by Owocki & Rybicki 
(1986) who found for a pure absorption LDW the sig- 
nal speed to be the sound speed and not the much 
faster Abbott speed. As noted already by these au- 
thors, this should be a consequence of assuming pure 
line absorption, which does not allow for any radia- 
tively modified, inward wave mode. Meanwhile there 
is ample evidence for Abbott waves in time-dependent 
wind simulations (Owocki & Puis 1999). 



Shallow solutions fail to reach infinity because 
they cannot perform the required spherical expan- 
sion work, implying that the flow starts to decelerate. 
Since this usually occurs very far out in the wind, the 
local wind speed is much larger than the local escape 
speed, and the wind escapes to infinity. Thus, a sim- 
ple generalization of the CAK model allowing for flow 
deceleration renders shallow solutions globally admis- 
sible. This raises a fundamental question of why the 
wind would adopt the critical solution at all, and at- 
tain the critical mass loss rate and velocity law, as 
proposed by CAK. 

In this Letter we analyze a physical mechanism 
which drives shallow solutions towards the critical 
one, and discuss under what conditions this evolu- 
tions does not terminate at the CAK solution, but 
continues into the realm of overloaded solutions. We 
find that simulations so far were affected by numeri- 
cal runaway towards the critical solution, by not ac- 
counting for Abbott waves in the Courant time step. 



2. Abbott waves 

Abbott waves are readily derived by bringing 
the wind equations into characteristic form. We con- 
sider a 1-D planar wind of velocity v(z,t) and density 
p(z,t), assuming zero sound speed. The continuity 
and Euler equation are, 
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Here, g(z) and F(z) are gravity and radiative flux, 
respectively. The CAK line force is given by g\ = 
C F (z) (v' / p) a (with v' = dv/dz), with constant C 
and exponent < a < 1. The unique, stationary 
CAK wind, v c (z),p c (z), is found by requiring a crit- 
ical point at some z c . The number of solutions for 
vv'(z) changes from 2 to 1 at z c (which is a saddle 
point), hence dE/d{vv') c = holds. Writing C in 
terms of critical point quantities, the Euler equation 
becomes, 
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Note that for stationary planar winds, pv is con- 
stant. If, in addition, g and F are taken constant 
with height, and p c v c v' / p is replaced by vv' /rh, with 
normalized mass loss rate rh = pv/p c v c , one finds 
that E does no longer depend explicitly on z for sta- 
tionary solutions. Hence, vv 1 is independent of z, 
too. This implies that z c is ill-defined, and every 
point of the CAK solution is a critical point. CAK 
removed this degeneracy by introducing gas pressure 
terms. Here we take a different approach and assume 
g = z/(\ + z 2 ). A situation with roughly constant ra- 
diative flux and gravity showing a maximum at finite 
height could be encountered above isothermal disks 
around compact objects (cf. Feldmeier & Shlosman 
1999). The critical point is determined by the regu- 
larity condition, dE/dz c = 0, hence z c = 1 and the 
critical point coincides with the gravity maximum. 
For simplicity also we chose a = 1/2 from now on, 
which is reasonably close to realistic values a < 2/3 
(Puis et al. 1999). None of our results should depend 
qualitatively on the assumptions made so far. The 
Euler equation is 



dv dv . . , dv/dz 

— +V—+ g(z) - J —!— = 0, (4) 

where g c = g(z c ). The stationary solutions for wind 
acceleration are given by 
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where plus and minus signs refer to steep and shallow 
solutions, respectively. For rh < 1, shallow and steep 
solutions are globally, i.e., everywhere, defined. For 
rh > 1, solutions are called overloaded, and become 
imaginary in a neighborhood of the gravity maxi- 
mum. These winds carry too large mass loss rates 
and eventually stagnate. 

Next we put the Euler equation into quasi-linear 
form, which does not mean to linearize it. Differenti- 
ating E with respect to z (Courant & Hilbert 1962; 



Abbott 1980) and introducing / = dv/dz, cqs. (0, |) 
become, 
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with inward Abbott speed in the rest frame, va = 
— y g c v/rhv' . In the WKB approximation, individual 
spatial and temporal variations are much larger than 
the inhomogeneous term g' / p in eq. (0), and the latter 
can be neglected. Consequently, v' j p is a Riemann 
invariant propagating at characteristic speed v + v a- 
Perturbations of v' / p correspond to the amplitude of 
a wave propagating at phase speed v + va- Note that 
v' j p is proportional to the Sobolev line optical depth, 
indicating that this wave is a true radiative mode. 

The second characteristic is determined by the 
continuity equation (^). In the advection operator 
in square brackets, v has to be read as v + in the 
zero-sound speed limit. This outwards propagating 
invariant corresponds to a sound wave, with ampli- 
tude p scaling with gas pressure. 

At the critical point, rh = 1 and vv'(z c ) = g(z c ) 
after eq. (|^), hence va c = —v c (where we introduced 
vac = va{z c ))- Abbott waves stagnate at the critical 
point, in analogy with sound waves at the sonic point. 
For shallow solutions, rh < 1 and vv' < v c v' c from 
eq. (H), hence v + va < 0. Shallow LDW solutions are 
therefore the subcritical analog to solar wind breezes. 

Because in the rest frame, the inward Abbott 
mode can propagate at larger absolute speeds than 
the outward sound mode, Abbott waves can deter- 
mine the Courant time step in time-explicit hydrody- 
namic simulations. Violating the Courant step results 
in numerical instability. Despite this fact, Abbott 
waves along shallow solutions were never considered 
in the literature. 



3. Wind convergence towards the critical 
solution 

We turn our attention to a physical mechanism 
which can drive LDWs away from shallow solutions, 
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and towards the critical one. Starting from an arbi- 
trary shallow solution as initial condition, we explic- 
itly introduce perturbations at some fixed location 
in the wind and study their evolution. In order to 
keep unperturbed shallow solutions stable in numeri- 
cal simulations, we fix one outer boundary condition, 
according to inward propagating Abbott waves. Ei- 
ther a constant mass loss rate at the outer bound- 
ary or non-reflecting boundary conditions (Hedstrom 
1977) serve this aim. At the inner, subcritical bound- 
ary, we also fix one boundary condition, according to 
incoming sound waves. Non-reflecting boundary con- 
ditions and p = const give similar results. 

Wind convergence towards the critical solution 
is then triggered by negative flow velocity gradients.^ 
Allowing for v' < turns the inward Abbott mode 
of phase speed v + v\ < in the rest frame into an 
outwards propagating mode. This is readily seen for 
a line force which is zero for negative v', i.e., when 
all photons are absorbed at a resonance location be- 
tween the photosphere and the wind point. The Eu- 
ler equation simplifies to that for an ordinary gas, 
with characteristic speed v — > in the zero sound 

speed limit. At the other extreme, for a purely local 

> 

line force where the unattcnuatcd stellar or disk ra- 
diation field reaches the wind point, g\ cx yjVj. Here 
the Abbott phase speed is found to be v + v\, with 
va = + \J —g c v /mv' for v' < 0. 

Consider then a sawtooth-like velocity perturba- 
tion (sinusoidal perturbation lead to similar results). 
Slopes v' > propagate inwards, slopes v' < propa- 
gate outwards. Hence, as a kinematical consequence, 
a sawtooth which is initially symmetric with respect 
to the underlying stationary velocity law evolves to- 
wards larger velocities. This is demonstrated in Fig- 
ure 1 where, in course of time, a periodic sawtooth 
perturbation is introduced at z = 2. The line force 
is assumed to be cx -\/|t>'|, and the initial shallow 
solution has rn — 0.8. The figure shows 2\ pertur- 
bation cycles. For upward pointing kinks the slopes 
propagate apart and a flat velocity law develops be- 
tween them. At each time step dt, a new increment 
dv = AdtSv/T (Sv and T being the amplitude and 
period of the sawtooth) is added at z = 2, hence the 
flattening velocity law does not show up in region A 



of Figure 1. Overall, the wind speed at the perturba- 
tion site evolves towards larger values during these 
phases. On the other hand, for downward point- 
ing kinks of the sawtooth, — Sv, the two approach- 
ing slopes merge, and the wind speed evolves back 
towards its unperturbed value after each decrement 
— dv — —AdtSv/T. The wind velocity hardly evolves 
during these phases, cf. region B of Figure 1. Over 
a full perturbation cycle, the wind speed clearly in- 
creases. 
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Fig. 1. — Upper panel: Evolution of a shallow wind 
during 2\ periods of a sawtooth perturbation of 
amplitude 10% being introduced at z = 2. Re- 
gions 'A' and 'B' correspond to phases where upward 
and downward pointing kinks, respectively, are intro- 
duced into the flow. Lower panel: stable Abbott wave 
excitation in the critical wind, at 35% perturbation 
amplitude. 

Essentially any perturbation which introduces 
negative v' will accelerate the wind. The amplitude 
of the perturbation is rather irrelevant since, with de- 
creasing perturbation wavelength, negative v' occur 
at ever smaller amplitudes. However, in more real- 
istic winds, dissipative effects may smear out short- 
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scale perturbations before they can grow. Details of 
the physical mechanism will be discussed elsewhere. 

//the perturbation lies downstream from the crit- 
ical point, the wind converges to the critical solution. 
Namely, as soon as the perturbation site comes to lie 
on the supercritical part of the CAK solution dur- 
ing its evolution, positive velocity slopes propagate 
outwards, and combine with negative slopes to a full 
wave train. No information is propagated upstream. 
This unconditional stability of the outer CAK solu- 
tion is shown in the lower panel of Fig. 1. 

4. Wind convergence towards overloaded 
solutions 

Wind runaway towards larger speeds as caused 
by perturbations introduced upstream from the crit- 
ical point does not terminate at the critical CAK 
solution. For low-lying perturbations, communica- 
tion with the wind base is still possible once the 
subcritical branch of the CAK solution is reached. 
The wind gets further accelerated into the domain 
of mass-overloaded solutions (where vv' > v c v' c and 
hence v > v c for z < z c according to eq. , until 
a generalized critical point develops, which prevents 
inward propagation of Abbott waves and adjustment 
of the mass loss rate. Such generalized critical points 
are given by 'termination' points, z t , of overloaded 
solutions, where the velocity becomes imaginary. At 
Zt, the number of real solutions vv'(r) changes from 2 
(shallow and steep) to 0. Hence, termination points 
arc defined by the same condition as the CAK criti- 
cal point (at which the number of solutions changes 
from 2 via 1 to 2), dE/d(vv')t — 0. From the sta- 
tionary version of eq. (^|), v^t — —ih, hence Abbott 
waves stagnate at termination points, and the latter 
are generalized critical points. 

The fact that perturbations with negative v' ac- 
celerate the wind either to the critical or an over- 
loaded state can be casted into black-hole conjecture 
(Penrose 1965): a LDW avoids a 'naked' base, and 
encloses it with a critical surface. 

Since to each z t there corresponds a unique, su- 



percritical mass loss rate, the latter is determined by 
the perturbation location alone. Using VAt — —vt, 
one finds rht = g c /gt > 1 for a planar wind with 
constant radiative flux. 

At a termination point, vv' jumps to the deceler- 
ating branch, vv' < 0. Beyond a well-defined location 
above the gravity maximum, the super-CAK mass 
loss rate can again be lifted by the line force, and 
vv' jumps back to the accelerating branch. Hence, 
two stationary kinks occur in the velocity law. Fig- 
ure 2 shows a hydrodynamic simulation of the evolu- 
tion towards an overloaded solution. Sawtooth-type 
velocity perturbations were introduced at z = 0.8. 
Correspondingly, m — 1.025 for the overloaded solu- 
tion, using g = z/(l + z 2 ). 




Fig. 2. — Wind evolution towards a stationary, over- 
loaded solution showing an extended decelerating re- 
gion. A periodic sawtooth perturbation is introduced 
into a shallow solution at z = 0.8, upstream from the 
critical point at z c = 1. 

Future work has to clarify whether LDWs show 
deep-seated perturbations. It seems however unlikely 
that they would occur at a unique location. Hence, 
overloaded winds should be non-stationary and show 
a range of supercritical mass loss rates. 

More fundamentally, time-dependent overloaded 
solutions occur already for single, unique perturba- 
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tion sites, once the latter lie below a certain height. 
For the present wind model, this is at z ~ 0.66. 
The overloading is then so severe and the decelerat- 
ing region so broad that negative wind speeds result 
(cf. Poe, Owocki, & Castor 1990). The corresponding 
mass loss rates are still only a few percent larger than 
the CAK value. The gas which falls back towards the 
photosphere collides with outflowing gas, and a time- 
dependent situation develops. Within each perturba- 
tion period, a shock forms in the velocity law, sup- 
plemented by a dense shell. These shocks and shells 
propagate outwards (Feldmeicr & Shlosman 2000). 

Although strong perturbations introducing neg- 
ative velocity gradients can appear already in O star 
winds, accretion disk winds are the prime suspects. 
The reason for this is that accretion processes and 
their radiation fields in cataclysmic variables and galac- 



verges to a critical or overloaded solution depends en- 
tirely on the location of perturbations, and not, e.g., 
on boundary conditions at the wind base. 

If Abbott waves are not accounted for in the 
Courant time step of hydrodynamic simulations, we 
find that numerical runaway can drive the solution 
towards the critical CAK wind. A detailed discus- 
sion of this will be given elsewhere. 

Future work has to clarify whether and where 
perturbations causing local flow deceleration, dv/dz < 
0, can occur in LDWs. Overloaded winds may be de- 
tected observationally. While their mass loss rates 
should still be close to CAK values, broad regions 
of decelerating flow could be identified in P Cygni 
line profiles. Furthermore, shocks occurring in over- 
loaded solutions with infalling gas may contribute 



to the X-ray emission from LDWs, besides shocks 
tic nuclei arc intrinsically variable on a range of timescalcs , ,. , . . .... /T ^ ^^ 

. ... „ „ . \ ....... „ „„. from the line-driven instability (Lucy 1982; Owocki 

et al. 1988). Note that the present wind runaway 



(Frank, King, & Raine 1992), and that disk LDWs 
are driven by a combination of uncorrected, local 
and central radiation fluxes. 



5. Summary 

We find that shallow solutions to line-driven winds 
are subcritical with respect to Abbott waves (sub- 
abbottic). These waves cause shallow solutions to 
evolve towards larger speeds and mass loss rates be- 
cause of the asymmetry of the line force with re- 
gard to positive and negative velocity gradients and 
because perturbations with opposite signs of dv/dz 
propagate in opposite directions. Steep velocity slopes 
propagate towards the wind base, steepen the inner 
wind and lift it to higher mass loss rates. In the pres- 
ence of enduring wind perturbations, this proceeds 
until a critical point forms and Abbott waves can no 
longer penetrate inwards. 

The resulting solution does not necessarily cor- 
respond to the CAK wind. For perturbations which 
originate below the critical point, the developing Ab- 
bott wave barrier is found to be the termination point 
of a mass-overloaded solution. The velocity law ac- 
quires a kink at the termination point, where the 
wind starts to decelerate. Whether the wind con- 



occurs already in the lowest order Sobolev approxi- 
mation, and is therefore, unrelated to the line-driven 
instability which depends on velocity curvature terms 
(Feldmeier 1998). 
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